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SIDESLIP IN A VISCOUS COMPRESSIBLE GAS* 
By V. V. Struminsky 


On the "basis of an analysis of the solutions of the equations 
of Navier-Stokes, it is found that in the slip of a wing in an iso- 
thermal or adiabatic flow of a compressible gas, the aerodynamic 
coefficients c x ^> C yl* and m z i can accura ' te - 1 -Y determined 

from the corresponding aerodynamic coefficients of the wing moving 
without slip. In the general case the aerodynamic coefficients of 
a sideslipping wing can be determined by the theory, developed in 
the present paper, of the three-dimensional boundary layer in a 
compressible gas . 


I. In the flow about a rectangular wing of infinite aspect 
ratio with sideslip, the velocity components u, v, and w, 
pressure p, density p, and temperature T in a system of coor- 
dinates attached to the wing (the z~axis parallel to the generator 
of the wing) will depend only on x and y. In this system of 
coordinates, the equations of Navier-Stokes and the boundary con- 
ditions can be written in the following form; 


ft du _ op xx ^ p xy 
p dt ' “5F" + 


Bpu , Spv _ 

ST 5T" 


dV _ ^Pyx 

P dt dx 
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p = pRT 
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where u = v = 0 on the surface of the wing; u = lbs cos 3 cos a 
and v = Voc cos p sin a at infinity; 



where w = 0 on the surface of the wing; w = Vo> sin 3 at infinity, 
3 is the angle of slip of the wing, and a is the wing angle of 
attack , 

In the case of isothermal flow (T = constant), the system of 
equations (la) can he satisfied by 

u = Vqd cos 3 Uq (x,y, Re co3 p, Ma cos p, a) 
v = Vod cos 3 Vq (x,y, Re cos 3, Ma cos 3, a) 
p = V ^ cos 3 p 0 (x,y, Re cos 3, Ma cos 3, a 

P = Pa> "Pq ( x >y> Ee 003 P> Ma cos 3> a ) 

The nondimensional functions Uq (x,y,Re,Ma,a) , Vg(x,y, Re,Ma,a), 
Pg(x,y, Re, Ma) and pg are a solution of the corresponding two- 
dimensional problem. The velocity component w(x,y,p) can be 
determined from equations (1c) and (2). 

In the case of adiabatic flow (p. = \ - 0), the system of 
equations (la) and (lb) can be satisfied by 
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u = Vo» cos 0 Uq (x,y, Ma cos _0,a) 

v = Voo cos 0 t q (x,y, Ma cos 0,a) 

P = P» V a. 003 3 Po (*>y> Ma 003 3, a) 

P = P 0 (x,y, Ma cos 0,<x) 


(3) 


In this case the velocity component w in the entire flow is 
equal to w = Voo sin 0 . 


II. In order to determine the components along the coordinate 
ares of the resultant force and resultant moment, the following 
expressions are used: 


X 1 = x l°( v ® coa 0 ,a) 

Yp = Y 1 °(V ot> cos 0,a) > 

M zl = M zl 0 (V oo cos 0 ,a) 


(4) 


Z 1 ds ^ 

Mxl lid) yds 


M y] _ = xds 




(5) 


dw 3w 

where (d = cos nx + ^ cos ny. 

Referring the forces and moments to the flow velocity at Infinity, 
for Isothermal flow. 


c xl = cos 2 

p^l 0 

(Re 

COS 

3, 

Ma 

cos 

0,a) 1 


Cyl = COS 

BOyl 0 

(Re 

cos 

3, 

Ma 

cos 

3, a) 

> (6) 

m zl = COs2 

n-m 0 
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cos 

3, 

Ma 

cos 

3, a) 
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For adiabatic flow, 


c xl = 

p 

cos^ 

P c xl° 

(Ma 

cos 

3, a) ^ 

1 

o 

II 

p 

003^ 

pc y i° 

(Ma 

cos 

P,a) 

> (?) 

m zl = 

p 

COS^ 

^zl° 

(Ma 

cos 

3> a ) ^ 



For adiabatic flow moreover, c zl = m^ = m^ = 0. 

Expressions (6) and (7) are direct consequences of the equations 
of Navier-Stokes . 


III. At large Reynolds numbers and at a Prandtl number equal to 
unity, the following system of equations for the three-dimensional 
boundary layer in a compressible gas is obtained from equations (la), 
(lb), and (lc) : 


f Bu Bu\ Bu B ( Bu\ -v 

p + v ^yj = - 5£ + SF^Syj | 



Bpu , Bpv _ 

Sy~ ~ 


0 


p = pRT 


> 


J 


(8a) 



a a 

= 



(8b) 


a Bw 
■ 37 “5F 


(So) 


where u = v = w on the wing surface; u - U(x, 1 'J c cos (3) and 
w = W= T 0= sinp on the outer limit of the boundary layer. 


The term U(x,Y<» cos 3) is the velocity of the adiabatic flow 
on the outer limit of the boundary layer. The value of U(x,Vcc ) 
can be determined by the method of 3. A. Christianovich (reference l). 
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In present work, in addition to the simplifying assumption 
Pr = 1, it is assumed that there is no heat transfer between the 
gas and the wing (bT/by)^ =0. In this case the simple integral 

u 2 + v 2 

of equation (8b) may be employed: ICpT = - = constant = 


By setting 
relation T = Tq 



V^p 

u 2 + w 2 \ 

2i 0 J 


as the temperature of stagnation, the 
is obtained. By the law of Bernoulli 
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p = p 0 



k 

u 2 + tr^Y ' 1 

21 o ) 


From the equation of state, 
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k 


c 0 
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U 2 ± W 2 Y -1 

2 ki 1 


(‘ 


2 2\ 
u + w \ 

21 0 ) 


For what follows, the following dependence of the viscosity 
on the temperature is assumed: p. = ^(t/Tq) 11 , where 0.5 <n< 1.5. 

We introduce the stream function pu = by /by , pv = -cty /bx and the 
new independent variables £ and T) : c£ /ox = p/Pqj dq/dy = p/Pq- 


A similar transformation was first applied by A. A. Dorodnitsyn in 
his work (reference 2). By carrying out the transformation of the 
variables in equations (8a) and (8c), the following system of equations 
is finally obtained: 


dv s 2 y by d 2 v 
ar ~ as ~ = 
br\ c 


d£3q cf ^2 P 01 - a Q 


1 - a - (3 


" 3o 


T# + 
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(1 - a - 3 ) 
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dq 
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where a = u 2 /21q, a q = U 2 /21q, 0 = w 2 /21 q, and 0 q = W 2 /21q. 

IV. The method of Kaman-Pohlhausen is used to solve the system of 
equations (9). After simple transformations and integration of the 
first equation with respect to from 0 to & and the second 
equation with respect to q from 0 to y , 



In these equations, 5 is the thickness of the boundary layer for 
the velocity component u, and y is the thickness of the layer for 
the velocity component w. 


By substituting in (ll) and (10) the expressions 


u _ 12 + \ AA _ \(A Z 
U " 6 2 ^5/ 




6 - \ 



w 

W 




( 12 ) 
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where X = 


a system of two ordinary differential 


PO S 2 dU/d j 

“ a 0 " 

equations of the first order is obtained for the determination of 
the boundary-layer thicknesses 5 and y. 


It is seen from expression (12) that in a compressible gas the 
direction of the velocity within the boundary layer does not coin- 
cide with the direction on the outer limit of the boundary layer. 
The flow is as if it were twisted within the boundary layer by a 
certain angle which changes with the height of the boundary layer 
and along the wing chord. Wear the point of separation in the 
lower layers of the boundary layer the flow is along the span of 
the wing. 

Y. For a flat plate about which there is a flow with sideslip, 
the system of equations (9) will have the following exact solution: 


ijr = 2|d ( 


■vs 


-Vcb cos |3£t 0 (T ) 


w = Vcc sin 3 cp' q (t) 


(13) 


where T =\J cos p/q^q/2. The function cp q (t) . satisfies the 
ordinary differential equation 


_d_ 

dT 


jjl - a Q 2 qj f 0 2 (T)J n_1 X<f> " 0 (t) * = - 2 <P 0 (t)9" 0 (t) (l4) 

and the boundary conditions = 0<P* (0) = 0; q = t» cp ' (*> ) = 1. For 
P = 0, equation (13) goes over into the solution of Dorodnitsyn for 
the flat plate (reference 2). 


On the basis of the expressions here given, 


'xl 


4 <P" 0 (0) 
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where as follows from equation (14), depends on 

2 l/2(k-l)Ma 2 

0 ” 1 - l/2(k-l)Ma 2 


The value of the drag coefficient of a flat plate is determined 
from the expression 

n-1 

*P" 0 (°) f ic . i a'N 2 

c„ = C V 1 cos p + c„ n sin (3 = - 1 + — 5 — Ma coa P 

X Xi Z1 V Re c* 008 P V / 
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